An ab-initio investigation of magnetism in two-dimensional uranium systems 
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The orbital and spin magnetic moments, and the X-ray-magnetic circular-dichroism (XMCD) 
spectra at the Mi^ edges of the U atoms in a UAs/Co multilayer and in an a-U film are calculated 
within the framework of the density-functional theory in combination with the local-spin-density 
approximation (LSDA), the generalized-gradient approximation (GGA) and the LDA+(7 method. 
The antiparallel arrangement between the U and Co spin magnetic moments at the interface results 
in the vanishing of ferromagnetism for the case of very thin Co layers. The U moments decay rapidly 
with the distance from the film surface. The magnitude of the magnetic-dipole term (T z ), which 
appears in the spin XMCD sum rule, is small. The different exchange-correlation treatments do not 
yield qualitatively different results. 
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Both the subject of actinides and the issue of mag- 
netism in two dimensional systems, like multilayers or 
films, represent interesting topics in contemporary solid- 
state physics. The joint problems from the two fields 
are therefore even more challenging. Actinides and ac- 
tinide compounds exhibit very different magnetic behav- 
iors, such as Pauli paramagnetism, localized and itin- 
erant magnetism and heavy fermionsi^ It is expected 
that the properties would change when going from bulk 
to two-dimensional materials. For example, Plasket et 
ali produced UAs/Co multilayers in an attempt to in- 
crease the Curie temperature of the amorphous UAs/Co 
alloy due to the induced polarization of U atoms through 
the exchange coupling to the Co layers. The first results, 
which came from magneto-optical measurements, were 
promising, but they were not reproduced with polarized- 
neutron-reflectivity experiments^ because the main con- 
tribution to the measured signal arises from the Co 
atoms. As an alternative approach, Kernavanois et alA 
performed a measurement of the X-ray-magnetic-circular 
dichroism (XMCD), which is element specific^. The same 
technique may be suitable for detecting the finite mag- 
netic moment at the surface of uranium that was pre- 
dicted theoretically^. 



The XMCD spectroscopy is based on the difference 
in the absorption coefficients /i + (e), At - ( e ) and (i°(e) for 
circularly- right, circularly-left and linearly polarized X- 
rays. The 5/ contributions to the magnetic orbital mi = 
~^b{Iz) and spin m s = — /is(cr z ) moments are related 
to the absorption coefficients at the M^g edges, via the 



sum rules&S: 
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and the total ab- 
1 The symbols 



with the XMCD signal fi c = /i + — fi~ 
sorption coefficient /i t = /i + + /i~ + 
(A t c,t)M 4i5 in Eqs. (3)-(5) denote that the integration is 
to be performed over the energy integral corresponding 
to the Mi t s edges. There are some problems connected 
with the application of the sum rules on experimental 
data. The exact number of holes ATj, in the / band is not 
directly measurable, nor is the expectation value (T z ) of 
the magnetic dipolar operator: 



£ = i[<7-3*(*.<r)] s 



(6) 



where cr denotes the vector of the Pauli matrices. The 
integrations in Eqs. (3)- (5) have to be performed over 
the energy interval which corresponds to the electron 
3c? — > 5/ transition. It is important to note that the 
limits of this energy interval are not uniquely defined in 
the experiment and that it is assumed that the electron 
3c? — > 5/ transition is the only relevant contribution to 
the absorption signal in this energy range. Above all, the 
sum rules are derived for the atomic case. In spite of all 
these limitations, the XMCD analysis, particularly at the 
^2,3 edges, has been proven to be a very successful tool 
for investigations of the magnetism in transition-metal 
systems of various dimensionalities. However, there is 
less research about the validity of the sum rules for other 



materials, for example, actinides, especially in systems 
other than bulk. Hence, the subject of the present paper 
is a comparison of the magnetic moments calculated di- 
rectly from the electronic structure with the moments ob- 
tained from the theoretical XMCD spectra by using Eqs. 
()1I2[) . The theoretical approach has several advantages 
over an experiment since some of the above-mentioned 
uncertainties, particularly ones related to the number of 
holes, the (T z ) term and the integration limits, can be 
reduced to a minimum. This enables us to investigate 
the validity of the sum rules without the influence 
of other effects present in the experiment, which are hard 
to control. The number of holes Ah is well defined. The 
integrations in Eqs. (3)-(5) are carried out over the in- 
terval between the Fermi energy E-p and the upper limit 
Eq , which is related to the and to the /-resolved den- 
sity of states n,f(e) as = J^ c n/(e) de. The (T z ) term 
is simply calculated as the expectation value of (6). 

The calculations were performed within the frame- 
work of the density-functional theory by applying 
the Wien97 codei^, which adopts the full-potential 
linearized-augmented-plane-wave (FLAPW) method 11 . 
The Brillouin-zone (BZ) integrations were carried out 
with the modified tetrahedron method^ by using 600 
and 550 k points for the UAs/Co multilayer and the U 
film, respectively. The plane-wave cut-off parameter was 
set to 11.4Ry for both systems. In order to investigate the 
influence of correlation effects, the exchange-correlation 
contribution to the effective potential was calculated by 
using different schemes, namely the local-spin-density ap- 
proximation (LSDA)i£, the generalized-gradient approx- 
imation (GGA)i4 and the LDA+J7 methodic. The pa- 
rameters U and J, which appear in the LDA+t/ scheme, 
were set to U — 2 eV and J — 0.5 eV^£. The contribu- 
tion of the spin-orbit coupling (SOC) was calculated in 
the second-variational schem o 18 ' 19 , since it was demon- 
strated that a more rigorous treatment, for example, the 
first- variational scheme^, did not yield substantially dif- 
ferent results for uranium, although the magnitude of the 
5/ uranium SOC interaction is comparable to the corre- 
sponding bandwidth. The absorption coefficients /z + (e), 
/i~(e) and /x°(e) as a function of the photon energy e were 
calculated using Fermi's golden rule in a nonrelativistic 
dipole approximation that is based on the evaluation of 
the matrix element for the operator p-e, with e denoting 
the polarization vector of the lighljSiiSSiSi. However, note 
that the current implementation of the density functional 
theory cannot simulate the excited states which would 
result from a hole in the 3d state. 

The UAs/Co multilayer was modelled with a supercell, 
based on a face-centered-cubic (fee) stacking sequence 
along the (001) direction as presented in Fig. 1. The ex- 
perimental lattice parameter a = 0.355 nm for fee Co was 
used. It turns out that the calculations performed on su- 
percells with fewer Co layers do not result in a magnetic 
solution, in agreement with experiments at 300 K (Ref. 
S) , where no magnetic ordering was observed for the mul- 
tilayers with the Co layer thinner than 1 nm. However, 



0.084, 1.62 



0.078, -0.16 1 



0.015,0.39 




0.075, 1.49 



U 



0.010, 0.40 



O As 



Co 



FIG. 1: A part of the As/Co- multilayer supercell, which is 
based on the fee Co unit cell with the lattice parameter a. The 
numbers (l z ), (a z ) represent the spin and orbital moments at 
particular U and Co sites obtained by applying the LDA+(7 
method for the 5/ U states. 



the present theoretical results do not prove the thesis 
that the absence of magnetism is due to the intcrdiffu- 
sion of U and As atoms into the Co layer. Furthermore, 
in the case of larger supercells, with more UAs layers, 
a finite magnetic moment at the U sites was found just 
at the interface, while in experiments^ the magnitude of 
the moment increased with an increasing UAs thickness. 
The calculated spin (a z ) and orbital (l z ) magnetic mo- 
ments, obtained with the LDA+C7 method, are presented 
in Fig. 1 . As demonstrated in Table 1 for the case of the 
U 5/ moments, the LSDA and GGA do not give qualita- 
tively different results, although the Coulomb repulsion 
U, in terms of the LDA+t7 method, enhances the or- 
bital moments by about a factor of 2 (3) relative to the 
LSDA (GGA) values. In contrast to the low-temperature 
experiments^ on the multilayers with amorphous UAs bi- 
layers, carried out in the presence of high magnetic fields, 
the magnitude of the orbital moments arc smaller than 
the magnitude of the spin moments. The magnitude of 
the total magnetic moment \{l z ) + (cr z )\ ~ 0.1 is much 
smaller than the corresponding experimental values 5 be- 
tween 0.4 and 1.1. The discrepancy might be partly 
due to the difference in the ratio of atomic concentra- 
tions, [U]/[As], between the calculation ([U]/[As]=l) and 
the experiment ([U]/[As] = 1.5) since the results of the 
SQUID measurements^ exhibited a pronounced compo- 
sition dependence. At the interface there is an antipar- 
allel alignment between the calculated U and Co spin 
moments. The Co moments in the interface layer are 
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FIG. 2: The 5/ orbital (l z ) (circles) and the spin (a z ) 
(squares) magnetic moments as a function of the layer num- 
ber for U film, calculated with the LSDA (upper graph) , GGA 
(middle graph) and LDA+U method (lower graph). 



larger than the magnitude of the orbital moment (l z ). 
The enhancement of the magnitude of (l z ) due to the 
Coulomb repulsion is smaller than in the case of the 
UAs/Co multilayer, whereas the LDA+[7 spin magnetic 
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drastically reduced relative to the moments in other lay- 
ers or in the bulk, which is consistent with the vanishing 
of magnetism in the case of thin Co layers, as mentioned 
earlier. 

A comparison between the U 5/ moments calculated 
directly from the electronic structure and the moments 
obtained from the theoretical XMCD spectra is given in 
Table 1. The orbital sum rule holds well for all three 
types of the exchange-correlation potential. The valid- 
ity of the spin sum rule J2J) is worse, since the resulting 
moments underestimate the directly calculated values by 
- 20% (LDA+J7) to - 30% (GGA). The influence of 
the magnetic dipole (T z ) term is below the deviation of 
the sum rule (a z ) from the directly calculated value al- 
though, surprisingly, the deviation is reduced when the 
(T z ) is omitted from Eq. J5J) in the case of the LSDA and 
LDA+Z7. 

The supercell for modelling the a-U (001) film was 
based on the orthorhombic structure and it contained 6 
uranium and 4 vacuum layers (for details, see, for exam- 
ple, Ref. 0). As found previously^, the density- functional 
theory predicts a magnetic ordering on the U surface. 
Fig. 2 demonstrates that the calculated magnetic mo- 
ments rapidly decrease with an increasing distance from 
the surface. However, the sign of the moments alternates, 
so that the moments in subsequent layers are aligned an- 
tiparallel to each other. In the surface and in the next 
layer, where the moments are possibly measurable, the 
magnitude of the spin moment (a z ) is about two times 



FIG. 3: The calculated absorption spectra for two polariza- 
tions (solid lines) and the XMCD signal (dashed lines) for 
the surface (upper graph), the second (middle graph) and 
the third (lower graph) layer of U film obtained by using the 
LDA+C/ method. The M5 edge is shown on the left and the 
M4 edge on the right panel. 



moment is almost two times larger than the correspond- 
ing LSDA value. The calculated XMCD signal is pro- 
nounced just for the surface layer, as shown in Fig. 3. 
Hence, only the surface magnetic moments are given in 
Table 2. Both the orbital and the spin sum rule under- 
estimate the directly calculated magnetic moments. The 
deviation is below 20%. The influence of the (T z ) term is 
even more subtle than in the case of the UAs/Co multi- 
layer, although the validity of the spin sum rule is worse 
if this term is set to zero. 

In conclusion, a theoretical test based on an ab-initio 
calculation proved the XMCD method to be a suitable 
tool for investigating the magnetism in two-dimensional 
U systems, that arises from the ordered 5/ magnetic mo- 
ments. However, the results of the calculations on the 
UAs/Co multilayer suggest the ferromagnetism confined 
at the interface, and hence to some extent disagree with 
the experimental data and findings. If ferromagnetism re- 
ally exists on the surface of uranium, its presence should 
be readily proven by means of the XMCD measurements. 

We are grateful to N. Binggcli for the reading of 
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* Electronic address: matej.komelj@ijs.si 



M. S. S. Brooks, B. Johansson, and H. L. Skriver, in Hand- 



4 





{Q 


(l z ) from Eq. JTJ 




(a z ) from Eq. © 


(a z ) from Eq. ©, (T z ) = 


LSDA 


0.033 


0.034 


-0.138 


-0.105 


-0.109 


GGA 


0.025 


0.027 


-0.169 


-0.121 


-0.104 


LDA+C/ 


0.078 


0.068 


-0.162 


-0.128 


-0.138 


TABLE I: A comparison between the 5/ contribution to the orbital (l z ) and the spin (a z ) magnetic moments at the U sites in 
the UAs/Co multilayer, calculated directly from the electronic structure, and obtained from the theoretical XMCD spectra by 
using the sum rules 11I2H. 
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(a z ) from Eq. 


(a z ) from Eq. ©, (T z ) = 


LSDA 


-0.444 


-0.364 


0.869 


0.722 


0.714 


GGA 


-0.576 


-0.479 


1.193 


0.981 


0.956 


LDA+U 


-0.655 


-0.546 


1.456 


1.194 


1.162 



TABLE II: A comparison between the 5/ contribution to the orbital (l z ) and the spin (a z ) U magnetic moments at the surface 
of the U film calculated directly from the electronic structure, and obtained from the theoretical XMCD spectra by using the 
sum rules 11121 . 
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